Studying Superfluid Transition of a Dilute Bose Gas 
by Conserving Approximations 
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We consider the Bose-Einstein transition of homogeneous weakly interacting spin-0 particles based 
on the normal-state ^-derivable approximation. Self-consistent calculations of Green's function 
and the chemical potential with several approximate $'s are performed numerically as a function 
of temperature near Tc, which exhibit qualitatively different results. The ladder approximation 
apparently shows a continuous transition with the prefactor c = 2.94 for the transition-temperature 
shift ATc/T^ = can^^^ given in terms of the scattering length a and density n. In contrast, the 
second-order, particle-hole, and fluctuation-exchange approximations yield a first-order transition. 
The fact that some standard <I>'s predict a first-order transition challenges us to clarify whether or 
not the transition is really continuous. 



I. INTRODUCTION 

The Bose-Einstcin condensation (BEC) of homoge- 
neous weakly interacting Bose gases has attracted much 
attention over a decadeF^r— As shown by Baym et al^^ 
this topic is profound enough to require treatments be- 
yond the simple perturbation expansion. To be specific, 
they confirmed that the transition temperature Tc starts 
to increase linearly with the s-wave scattering length a 
as 



,1/3 



(1) 



where n is the density, and presented analytic esti- 
mates for the prefactor c using various approximations 
for Green's function. Subsequently, a couple of Monte 
Carlo simulations on finite lattices obtained a widely ac- 
cepted value c « 1.3i^ 

However, these studies as well as othersir— focused 
mostly on the critical density by implicitly assuming a 
continuous transition, thereby leaving behind an impor- 
tant question of how the system approaches the critical 
point as a function of temperature. 

We will consider the issue based on the conserving $- 
derivable approximationji^^— This systematic approxi- 
mation scheme has several remarkable advantages,^ as 
may be realized by the fact that the Bardeen-Cooper- 
Scheriffer theory of superconductivity belongs to it as a 
lowest-order approximatioujSS and has been used exten- 
sively to clarify anomalous properties of high-Tc cuprate 
superconductors.'^*' Thus, the method will help us to see 
the critical region of BEC more closely. Indeed, the con- 
tents here may be regarded as an extension of those with 
self-consistent approximations by Baym et al^i^ just on 
T = Tc to (i) incorporate temperature dependences of 
T ^ Tc and (ii) consider more approximations systemat- 
ically. We will thereby find that the self-consistent one 
bubble approximation they considered yields a first-order 
transition contrary to their assumption. The main results 
are summarized in S3. 2 below. 



II. FORMULATION 

A. Hamiltonian and Green's function 

We will consider identical homogeneous bosons with 
spin 0, mass to, and density n interacting via a weak 
contact potential U5{ri — r2). To study this system near 
Tc, we adopt the units 
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where C(3/2) = 2.612 •• • is the Riemann zeta function 
and /cb denotes the Boltzmann constant. Thus, the 
critical temperature of the ideal Bose gas,^'^'^^ T° = 
(27r;i^/fcBTO)[n/C(3/2)]2/3, is set equal to 1, and the ki- 
netic energy is expressed in terms of the momentum p 
simply as Cp = ■ 



The Hamiltonian is given by 
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where /i is the chemical potential and c^, and Cp are the 
creation and annihilation operators, respectively. Ultra- 
violet divergences inherent in the continuum model are 
cured here by introducing a momentum cutoff S> 1. 
However, our final results will be free from pc, as seen 
below. It is standard in the low-density limit to remove 
U in favor of the s-wave scattering length a. They are 
connected in the conventional units by^^ 
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d^p 6{pc-p) 



with 9(x) the step function, which in the present units 
reads l/Svra ^ 1/U +Pc/47r^. We will focus on the limit 
a — and choose Pc so that 1 pc 7T/2a is satisfied. 
Thus, we can set 



U = Stto 
to an excellent approximation. 
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FIG. 1: Diagrammatic representation of symmetrized vertex 
r'''-'(pi,p2, (7), which is equal to 2(7 for the contact interaction 
U5(r\ — r2) used here. 



Let us introduce Green's function in the normal state 



by 



equal to 2/7 with no momentum and frequency depen- 
dence for the present contact interaction U5{ri — r2).— 
To begin with, ^iflex is given as a sum of four kinds 
of diagrams in Fig. [2] as 



^•flex = $1 + $2 



ph 



(8) 



where $1 = VUv? is the Hartree-Fock term. To express 
the other contributions analytically, let us introduce the 
functions 



.{q)^^.Y.G{p)G{p + q), 
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where p = {p, i£„) with e„ = 2mrT (n = 0, ±1, ±2, • • • ) 
a boson Matsubara frequency. The self-energy S(p) is 
given exactly by 



1 (5$ 

'T5G{py 



(5b) 



where the functional $ = *i'[G'] is defined as the infinite 
sum of closed skeleton diagrams in the simple perturba- 
tion expansion for the thermodynamic potential with the 
replacement Gq Gii^iiSiSi The chemical potential ji is 
connected with the particle density n by 



n = -^^G(p)e-"°+, 



(6) 



Xpp(9 ) ^T,Y1 G{p)G{^p + g ), 



V 



(9b) 



each of which corresponds to a pair of lines connect- 
ing adjacent vertices in Fig. [2](c) and (d), respec- 
tively. By following the Feynman rules for the pertur- 
bation expansion in terms of F^"^ = 2C/,^2, the numer- 
ical factor of Fig. [UJb) and those of the nth-order di- 
agrams {n = 3, 4, • • • ) in Fig. and (d) are easily 



found as C2 



-2Uf2'^T/e2\ ^ -C/2T/2, 



.(") 



-(-2J7)"22"-i(n - l)!T/4"n! (-2)"-i;7"T/n, and 

= -(-2C/)"2"(n-l)!T/4"n! = (-l)"-iC/"r/n. We 
thereby obtain analytic expressions for Fig.[Ub)-(d) as 



T 



T 



E[t^xph('f)]' = -^E[f^^pp(9)]'' (10^) 



with V the volume and 0+ an infinitesimal positive con- 
stant. 

If BEG is realized as a continuous transition, the tran- 
sition temperature will be determined by the condi- 
tion: 



T 

*ph I'^tl + 2C^Xph(g )] - 2C/xph(g ) 



Pi = S(0). 



(7) 



This relation, which is derived from the Hugenholtz-Pines 
relation in the condensed phas O ' ^^'^^ by setting the off- 
diagonal self-energy zero, naturally extends the condition 
/io(l) = for the transition temperature of the ideal gas^i 
to interacting cases. Indeed, it was used by Baym et al^^ 
to estimate the critical density Uc- 



B. FLEX and related approximations 

The <I>-derivable approximation consists of (i) retain- 
ing some finite terms or partial series in $[G] and cal- 
culating G and S self-consistently by eqs. ([5a|) and 
(jSbp . We will consider the fiuctuation-exchange (FLEX) 
approximation^- and those derivable from it by reducing 
terms. All of them are concisely expressible in terms of 
the symmetrized vertex T^'^\pi,p2,q) of Fig. [TJ which is 
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[Uxppiq) 



(lOb) 



(lOc) 



respectively. 

Let us substitute eq. dH) into $ of eq. ([Sb]) . We then 
obtain the self-energy as 

E(p) ^2Un-y Y.^V2{q) + Vpy,{q)]G{p + q) 
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with 
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V2{q) = -2U\pi,{q), 
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(a) Hartree-Fock 



(b) 2nd 






(c) particle-hole (d) particle-particle 

FIG. 2: Feynman diagrams for "I>flex 



[2t/Xph(g)]^ 
l + 2C/xph('?)' 



Vppiq}=2U 
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Equation ([5a|) with eq. (fTT|) forms a closed nonlinear 
equation for G{p) that may be solved numerically to clar- 
ify normal-state properties of T > Tc- 

Besides eq. (fTTj) . we will also consider the following 
self-energies: 



I]2(p1 = 2C/n - ^ ^ V2{q)G{p + q), 



V 



(13a) 



Eph(p) = 2C/n - ^ ^[1^2(9) + ^ph(«)]G(p + g), (13b) 



V 



T 



Spp(p) ^ 2C/n - - Y,[V2{q)G{p + q) 

q 

+Vpp{q)G{-p + q)], 



(13c) 



which are all derivable from eq. (jlip by reducing terms. 
Using these different approximations, we may check how 
well the present approach describes the BEC transition. 

Equations (|13a[) and (|13cP were used by Baym et al.—i^ 
as "self-consistent one bubble approximation" and "self- 
consistent ladder sum" for estimating Tc. They also 
considered the "self-consistent bubble sum" of setting 
^ph{q) — ^Vph{q) in eq. (jl3bj) . which amounts to neglect- 
ing the exchange processes altogether in the particle-hole 
series of Fig. ^c). Since the exchange processes are def- 
initely present, we will consider eq. (jl3bp instead of the 
self-consistent bubble sum. 



C. Dilute gas near Tc 

We now focus on the critical region of the weak- 
coupling limit a — 0, i.e., T > 1. Noting that this system 



is quantitatively close to the ideal gas, we first express 
the chemical potential as 



= /^o + S(0) + A^, 



(14) 



where /jq is the chemical potential of the ideal gas van- 
ishing quadratically for T — 1 as^ 



Co 



3C(3/2) 
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We also adopt the classical-field approximation by Baym 
et al^^ of retaining only the e„ = component in eq. 
([5a| : its validity will be confirmed shortly. We subse- 
quently perform a change of variables, T — ^ r, p — >■ p, 
A/i — > r, and S(p, i£„ = 0) — ?> cr(p), given explicitly by 



T = l + ^., 
Co 



p = U Tp, 



Mo(T) + AAi(r) = -(;7T)2r(r), 



E(p,ie„ = 0) = S(0) + (UrfaiP). 



(16a) 



(16b) 



(16c) 



(16d) 



Note that F — 1 and ct(p) for [/ — 7> as seen from 
eqs. (|15l) and The change of variables also removes 
the remaining source of the ultraviolet divergence, i.e., 
S(0), completely from the theory. Substituting eq. 
into eq. ((5a|) for e„ = and using eq. p6)) . we can write 
G(p, iSn = 0) as 



G(p,0) 



g(p) 



ff(p) 



p^ + ct(p) + F 



(17) 



We also adopt the classical-field approximation for eq. © 
and substitute eqs. (I16al) . ()16bl) . and (|T7)) into it. It then 
turns out that Xph(g,0) = Xpp(9:0) = x{q)/Ut with 



x{q) 



d^P 
(2^ 



5(p)5(|p + g|)- 



(18) 



Subsequently using eqs. PT|) and (jl6dp with eqs. (jl6al) . 
(|16b|) . and PT)) . we obtain an equation for the reduced 
self-energy a{p) as 



cr(p) 



d3g 

(27r)^ 



-«(g)[5(|P + 'Z|)-<?('Z)], (19) 



with 

D(g) 



X(g) , [2x(g)/T] 



[X{q)/r? 

l + 2x(g)/r ^ l + x('7)/7 



(20) 



As for eq. ^ for the chemical potential, we subtract 
its non-interacting correspondent from it, adopt the 
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classical-field approximation, and use eq. p7|) . It is 
thereby transformed into an equation for F = r(r) as 
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p2 + (t(p) + r p'^ + 1 



p'^dp - 



0. 



(21) 



Equations (fT9)) and ([2T|) form closed equations for a(jj) 
and r for a given r. Note that the transformation (fT6|) 
has removed completely from the self-consistent equa- 
tions. Besides, they are free from ultraviolet divergences. 

The self-energies of eq. can be transformed simi- 
larly, which turn out to have the kernels 



_ 2x{q) 1 - 2x(g)/T 

1 + 2x(g)/T 



(22a) 



(22b) 



(22c) 



in place of v{q) in eq. (|T9)) . Note that ^2(9) and Wpp(9) are 
negative, as seen from eq. (IT51) . whereas ^(g) and 'yph('z) 
change sign from negative to positive as r is decreased to- 
wards zero. These distinct behaviors of the kernels from 
different approximations will lead to contradictory pre- 
dictions on the BEC transition, as seen below. 

Finally, eq. (O for the continuous transition point is 
transformed with eqs. (HH) and (jl6cp into 



r(Te) = 0. 



(23) 



The corresponding critical temperature is easily obtained 
by eq. p6a|) with r = Tc. Using eqs. jH), (g]), and ^ as 
well as = 1 in the present units, we confirm that the 
transition-temperature shift starts linearly in a as eq. ([T|) 
with the prefactor 



3[C(3/2)]4/3 



Tr = 58.52t^ 



(24) 



A couple of comments are in order before closing the 
subsection. First, eq. tells us that diagrams from 
second through infinite orders in U contribute equiva- 
lently to cr(p). The validity of this statement is clearly 
not restricted to the FLEX approximation alone; it can 
be confirmed easily by applying eqs. (ITBl) and (IT71) for 
general nth order terms in the classical-field approxima- 
tion. Thus, we need to include infinite diagrams of $ 
to obtain an exact value of Tc in the self-consistent per- 
turbation approach, which is practically impossible. We 
may expect, however, that some approximations for $ 
enable us to obtain qualitatively correct results for the 
BEC transition. Second, e„ 7^ components in eq. ([5a|) 
are smaller than the e„ = one by {Ut)'^ in the critical 
region so that they are negligible, as seen easily by using 
the transformation of eq. (|16p . Thus, the classical- field 
approximation by Baym et al^f^ has also been justified 
by the present consideration. 



III. RESULTS 
A. Numerical procedure 

We explain how to solve eqs. (|19l) and (PT|) numerically 
to obtain the reduced self-energy a{p) and reduced chem- 
ical potential F as a function of the reduced temperature 
T. Let us introduce the non-interacting correspondent of 
eq. ((IHl) as 

d?p 



1 



(27r)3(|p + q/2|2 + l)(|p-q/2|2 + l) 
arctan(q/2) 



A'Kq 



(25) 



The second expression has been obtained by (i) per- 
forming angular integrations, (ii) subsequently expand- 
ing ln[(p± 9/2)2 _^ ^ ievTos of pq/ {p^ + / I), (in) 
making a change of variables as p = ((7^/4 + l)^/-^tan6' 
and carrying out the 9 integration, and (iv) comparing 
the resulting series with the Taylor expansion of arcsinx. 
Noting that x(g) X^^H?) foi' 9 — >■ 00, we reahze that 
x{q) — > l/8g for 9 ^ 1. Hence, it follows that eq. ([19]) 
for p ^ 1 is well described by 



.(0) 



dh 



1 



(27r)3 4r2g 
ln(p2 + 1) - 2 



1 



q^ + 1 \p + q\ 
- (2/p) arctanp 



(47rr)2 



(26) 



where the second expression has been obtained in the 
same manner as eq. (j25p . 

Equations (jl9p and (f2T|) have been solved iteratively, 
starting from the non-interacting Green's function and 
chemical potential in the integrands. Functions x(g) and 
cr(p) are calculated by using the integral expressions of 
x{q) ~ X^^\q) and a{j)) — CT^°^(p) whose integrands de- 
crease more quickly in the high-momentum region than 
those of eqs. (IT8|) and ([19]). Further, we make a change of 
variables p — sinh v? and q — sinli for the integrations 
to cover a wide momentum range up to a cutoff momen- 
tum pcut = 10^ ~ 10^. We have also stored x{q) and a{p) 
at equal intervals in terms of v and u. These values are 
used in the next step of iteration with interpolation. The 
region p -I- q of eqs. (fTS)) and (|T9]) are handled sepa- 
rately to incorporate more integration points in both the 
polar and radial integrations. The convergence has been 
checked by changing the number of integration points as 
well as pcut- 

The above procedure has been repeated for four differ- 
ent approximations with kernels (|20)) and ([22]) to check 
how reliable the predictions on the BEC transition by the 
present approach are. 



B. Results 

Figure [3] plots reduced chemical potential F cx S(0) — /z 
as a function of reduced temperature r cx T — 1 in the (a) 
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FIG. 3: Reduced chemical potential F as a function of re- 
duced temperature r in (a) the self-consistent second-order 
approximation and (b) self-consistent ladder approximation. 



self-consistent second-order approximation with kernel 
(|22ap and (b) self-consistent ladder approximation with 
kernel (|22cp . A continuous BEC transition corresponds 
to a monotonic decrease of T{t) towards where BEC is 
realized. Thus, Fig. El^b) from the ladder approximation 
apparently exhibits a continuous BEC transition around 
Tc = 0.0503, which translates to ATc/T° = 2Man^/^ by 
eqs. (P and ([M]). On the other hand, Fig.E^a) from the 
second-order approximation shows a clear sign of a first- 
order transition somewhere between 0.074 ^ r ^ 0.076 
where r(r) is multivalued. As for the FLEX approxima- 
tion, we have not even found a solution that approaches 
F = continuously; here F(t) has a minimum F,„iii 0.95 
around r = 0.14. This strange behavior is brought about 
by kernel (PO)) that changes sign from negative to pos- 



itive as T is reduced. The same statement holds for 
the particle- hole approximation of eq. (I22b|) . Since the 
present model is quite close to the ideal Bose gas, we may 
conclude that the BEC transition is a first-order transi- 
tion in both the FLEX and particle-hole approximations. 

With the diversity of predictions in the present ap- 
proach, we can hardly say anything definite about the 
nature of the BEC transition with a weak two-body in- 
teraction. However, the fact that some standard ^'s (i.e., 
$2, + and $FLEx) exhibit a first-order transition 
challenges us to clarify unambiguously whether or not 
the BEC transition is really continuous. 

As for the transition temperature shift, the value 
ATc/T° = 2.94an^/'^ in the present ladder approxima- 
tion departs from the value ATc/T°~2.5ari^/^ obtained 
by Baym et al. ^ with an ultraviolet cutoff A. How- 
ever, they also reported that their bubble-sum result 
ATc/T° ~ 1.6an^/'^ for the finite A is increased up to 
2.0an^/'^ by the extrapolation A — oo. Incorporating the 
same difference 0.4an^/'^ into their ladder-summation re- 
sult yields ATc/TO ~ (2.5 + OA)an^/^, which is in good 
agreement with 2.94an^/'^ in the present approach. If 
we use Tc ~ 0.076 from Fig. [2{a) determined by eq. 
([23)) . we obtain ISTc/T^ = AAan^/^ for the second-order 
approximation, which is also in good agreement with 
ATc/T^ ~ (3.8 + 0.4)ani/3 by Baym et al^ However, it 
should be noted once again that our second-order result 
clearly indicates a first-order transition contrary to their 
assumption. In their early study, Baym et al.^ also re- 
ported ATc/T° = 2.9ani/3 3,3 eq. (29), which apparently 
agrees with our ladder result mentioned above. How- 
ever, their estimate starts from the second-order approx- 
imation corresponding to our eq. (|13ap and subsequently 

adopts a trial form U{k) « k^/^kc^^/[l + {k/k^f^^] for 
U{k) = 2m[T,2{k,ien = 0) - ^] at T = Tc to interpo- 
late between the low- and high-momentum behaviors of 
efc + S](fc, iSn — 0) — with an intermediate parameter kc- 
Thus, the agreement is accidental and their prefactor 2.9 
should be replaced by 3.8(-|-0.4) in their later numerical 
study for the second-order approximation^ 
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